EQUIDISTRIBUTION OF THE FEKETE POINTS ON THE SPHERE 
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ABSTRACT. The Fekete points are the points that maximize a Vandermonde-type determinant 
that appears in the polynomial Lagrange interpolation formula. They are well suited points for 
interpolation formulas and numerical integration. We prove the asymptotic equidistribution of the 
Fekete points in the sphere. The way we proceed is by showing their connection with other array 
of points, the Marcinkiewicz-Zygmund arrays and the interpolating arrays, that have been studied 
recently. 



1. Introduction 

For any integer I > 0, let He be the space of spherical harmonics of degree I in E> d . For any 
integer L > we denote the space of spherical harmonics of degree not exceeding L by ITl 
These vector spaces have dimensions 

,. „ d + 2L fd + L-l\ Td 

Let {Qi, . . . , Q% L } De any basis in U L . The points Z L = {zli, • • • , zl,w l } maximizing the 
determinant 

|A(xi, . . . ,x VL )\ = \det(Qf(x j )) i j\ 

are called the Fekete points of degree L for E> d (this points are sometimes called extremal fun- 
damental systems of points as in [SW04]). They are not to be confused with the elliptic Fekete 
points which are a system of points that minimize the potential energy. The extremal fundamental 
system of points are better suited nodes for cubature formulas and for polynomial interpolation, 
see [SW04] and the references therein. 

The geometric properties of the distribution of the Fekete points on the sphere has been the 
subject of research, see for instance [Rei90], [LBW08] or [SW04]. One natural problem is the 
limiting distribution of the points as L —> oo. If we denote by hl = ^- Ylj $xj ^ nas been known 
long ago that for the elliptic Fekete points p, L converges vaguely to the uniform distribution of 
the sphere even for a wide class of potentials, see [HS04] for a very nice survey. For the Fekete 
points in compacts K C C this is a classical result. Much less is known in higher dimensions. On 
the paper [BB08] the authors have found the limiting distribution in the context of line bundles 
over complex manifolds. The techniques there are very different from ours, they rely on a careful 
study of the weighted transfinite diameter and its differentiability. We will rather emphasize the 
connection of the Fekete points with the Marcinkiewicz-Zygmund arrays and the interpolating 
arrays (see below for the definitions). As long as the density of these arrys is understood we 
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can obtain the equidistribution of the Fekete points. This is the case of the sphere where we can 
build on the work [Mar07], where the M-Z arrays and interpolating arrays are studied. The same 
approach is being pursued by R. Berman in line bundles over complex manifolds. 

1.1. Marcinkiewicz-Zygmund inequalities and interpolation. We consider arrays of points 
on the sphere S ' that determine the norm of the polynomials, and also arrays of points where we 
are free to interpolate arbitrary values by polynomials,!. e. interpolating arrays. More precisely, 
For any degree L we take m l points in E> d 

Z{L) = {z Lj E$ d :l<j< m L }, L>0, 

and assume that m L — > oo as L — > oo. This yields a triangular array of points Z = {Z(L)} L > 
in S d . 

Definition 1.1. Let Z = {Z(L)}l> be a triangular array with > n l for all L. We call Z 
an L p -Marcinkiewicz-Zygmund array, denoted by L P -MZ, if there exists a constant C p > such 
that for all L > and Q eUl, 

ifl < P < °o 5 and 

sup \Q(uj)\ < C sup \Q(z Lj )\, 

ui& d j=l,...,m L 

when p = oo. 

Then the L p — norm in S> d of a polynomial of degree L is comparable to the discrete version 
given by the weighted £ p — norm of its restriction to Z(L). 

Definition 1.2. Let Z = {Z(L)}l>o be a triangular array with tul < ^ifar all L. We say that 
Z is L p —interpolating, if for all arrays {cLj}L>o,i<j<m L of values such that 

sup — \c Lj \ p < oo, 

L>0 7TL , 

there exists a sequence of polynomials Ql E TIl uniformly bounded in LP such that Ql^zlj) = 
clj, 1 < 3 < m L - 

Roughly speaking in order to recover the If— norm of a polynomial of degree L from the 
evaluation at the points in Z(L) we need a sufficiently big number of points in Z(L). Thus, 
intuitively, the M-Z arrays must have high density. On the other hand, in an interpolating array it 
is possible to have a spherical harmonic of degree at most L attaining some prescribed values on 
Z{L). Intuitively this is possible only when Z(L) is sparse. 

In dimension one the roots of unity are simultaneously an interpolating and an M-Z array when 
1 < p < oo. On higher dimension the situation is more delicate. It has been proved, [Mar07, 
Theorem 1.7] that there are not arrays which are simultaneously L P -MZ and interpolating when 
d > 2 and p ^ 2 and most likely even when p = 2. We will prove that the Fekete points are 
a very reasonable substitute. If we perturb them slightly they are interpolating sequences and 
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a different perturbation makes them MZ-arrays. Thus, in a sense, they behave like the roots of 
unity in higher dimensions. Since the densities of the MZ-arrays and the interpolating arrays are 
well understood, see [Mar07, Theorem 1.6], then we will get some geometric information on the 
Fekete points. 

In the next section we will provide the connection of Fekete points and interpolating and MZ- 
arrays. On the last section we will draw some geometric/metric consequences. 

Acknowledgment This paper has its origins in a conversation of the second author with Robert 
Berman at the Mittag-Leffler Institute over the possibility of connecting Fekete points and sam- 
pling sequences. It is a pleasure to thank him for sharing his thoughts and to the Institute for the 
warm hospitality and great atmosphere. 



2. Fekete points, MZ-arrays, Interpolating arrays 

Theorem 2.1. Given e > let L £ = [(1 + e)L] and 

Z £ (L) = Z{L e ) = {z Le ,i,...,z Lej7rLe }, 

where Z(L) is the set of Fekete points of degree L, then Z £ = {Z £ (L)} L > is an LP-MZ array, 
for any 1 < p < oo. 

Proof. Assume that Z is a collection of Fekete points. It satisfies a nice separation property that 
is convenient to prove the first inequality of (1.1). 

Definition 2.2. A triangular array Z is uniformly separated if there is a positive number s > 
such that 

e 

L + l 

for all L > 0, where d(z, w) = arccos(2;, w). 



d(z Lj , z Lk ) > 1 — T , if j ^ k, 



Reimer, [Rei90] observed that the Fekete points are uniformly separated. More precisely, 

7T . 

— < mm d(z Li ,z L j) 

just by using Marcel Riesz result for trigonometric polynomials on great circles. Thus we know 
that 

■At \ ^ * ^ ^ 

mm d{zL e i,z Lej ) > Tpr - TXT' 
ZL £ -L+l 

and therefore the array Z £ is uniformly separated. This implies the following Plancherel-Polya 
type inequality for any 1 < p < oo 

i 7ri£ r 

— / \Q(z)\ p M*), foranyQen^ 

see [Mar07, Corollary 4.6]. 

The right hand side inequality in (1.1) is more delicate, we need an appropriate representation 
formula for the polynomials in terms of the values at the points. The most naive approach is to 
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start by the Lagrange interpolation formula. Let 

n I \ _ A(z L1 , ■ ■ ■ , Zl.z-I; Z, ZL,j+l, ■ ■ ■ , Zl,tt l ) 
"Li[ z ) — T~/ s 

A(z L1 , . . .,Z L ,ir L ) 

then sup ze g d \£Li(z)\ < 1 and the Lagrange interpolation operator defined in C(S d ) as 

A L (f)(z) = J2f(zL J )h J (z) 

satisfies 

IIA^IIoc^^H/IU. 

We want better control of the norms. So we need a slightly bigger set of points and a weighted 
representation formula. Let p be a polynomial in one variable of degree [Le] and such that 
p(l) = 1. Then given Q e TIl one has for a fixed z E E> d 

R(w) = Q(w)p((z,w)) en Le 

and therefore we obtain our weighted representation formula: 

We define the operator Q L from C 71 "^ — > n i2£ as 

Ql[v]{z) = X>iK<*.*W>Kw(*) V ^ G 
i=i 

We want to prove that 

(2) / |Q^M P eM3)<— 5>,f, 

with constants uniform in L which is the righthand sided of (1.1). We need to choose the weight 
p with care. We need a polynomial p that peaks at one point, has degree [eL] and decays fast 
far away from the picking point. For this purpose we will use powers of the Jacobi polynomials 
which are natural in this context because they are the reproducing kernels in ILl, see [Mar07]. 
The Jacobi polynomials pj?'W of degree L and index (a, (3) are the orthogonal polynomials on 
[—1, 1] with respect to the weight function (1 — x) a (l + x) 13 with a, (3 > — 1. We take the 
normalization 



L 

We can use the estimates in [Sze75, Section 7.34] to obtain, for any v G S d 

(3) / \P¥ /2 ' d/2 - 1 \(u,v))\*da(u)~l, VL>0. 
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We will use as auxiliary polynomial 

then p(l) ~ 1 and by the estimate (3) 

/ \p((z jZL J)\da(z)^L- d [ \P^f- 1 \(z,z Le M 2 d*(z)^L- d ^-Kj; 1 . 
Js d Js d 

Now as \£L e ,j(z)\ < 1 one has 

/ |QlHI^(^)<^|^| / \p((z,z Le MM*)S— EKI 

and also for any fixed z e § d 

\Ql[v}(z)\ < sup\vj\^2\p((z,z Le j))\ 
i j= i 

<Bup|«j| / n L \p((z,z Lstj ))\da(z) <sup\vj\. 

3 JS d j 

Then the result follows by the Riesz-Thorin interpolation theorem. □ 

The corresponding result for interpolation reads as follows: 

Theorem 2.3. Given e > let L_ e = [(1 — e)L] and let 

Z-e{L) = Z(L^ e ) = {z L _ s>1 , z L _ Sj7rL _ £ }, 

where Z{L) is the set of Fekete points of degree L, then the array Z_ e = {Z„ e (L)} L > is 
LP -interpolating, for any 1 < p < oo. 

Proof. Given an array of values {v L _ e j} j= [ E , we can define the polynomials in H L 

r l[v]{z) = J2 v L- e jP{{z,z L _ e>j ))£ L _ Ej {z), 
i=i 

and Rh[, z L. E j) = v L_ e j- This time the map Rl is from C* L - e — > Hl and the L p -estimates on the 
norm of Rl follow exactly as the estimates of Ql in the previous Theorem. □ 

3. Geometric properties of the Fekete points 

We will draw some geometric information on the Fekete points. For a given z G S d and 
< R < 1 we denote by B(z, R) the spherical cap B(z, R) = {w e § d ; d(z, w) < R}. We will 
prove that as L — > oo the number of Fekete points in B(z, R) gets closer to iTLa(B(z, R)) where 
a is the normalized Lebesgue measure on S d , i.e. a = a/a(E> d ). We need first information on 
the density of M-Z and interpolation arrays. 
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Definition 3.1. For Z a uniformly separated triangular array in E* d we define the upper and 
lower density of the array respectively as 

mm z&d #{Z{L) n B(z, a/L))/n L 



D (Z) = lim inf lim inf 

a— >oo L— >oo 

D + (Z) = lim sup lim sup 



a->oo l^oo a(B(z,a/L)) 

min zeS d #(Z(L) n B(z,a/L))/n L 



a. — >oo L — >oo 

a(B(z, a/L)) 

The main result in [Mar07, Theorem 1.6] is (with a slight different notation) 

Theorem 3.2. Let 1 < p < oo. Let Z be a uniformly separated array. If Z is an LP- 
Marcinkiewicz-Zygmund array then D~(Z) > 1. On the other hand if Z is an LP -interpolating 
array then D + (Z) < 1. 

Let Z{V) be the set of Fekete points of degree L. We know that for any e > the array 
Z £ = {Z £ (L)} l>o is L 2 -MZ, so if we unwind the definitions corresponding to the densities, we 
get that for any e > 0, there is a big a = a(e) such that for all L and z G S d 

(4) nL — L > ri-e) 

Similarly since Z_ £ is interpolating, whenever Z is a Fekete array, from the density condition 
we get that there is a big a = a(e) such that for all L and z E E> d 

<5) — m^m — - (1+£) - 

3.1. Vague convergence. Let us see how the inequalities (4) and (5) imply that the normalized 
counting measure converges vaguely to the Lebesgue measure. Indeed, defining 

1 nL 

we have 

(fJ-L * Xb(n,<x/l))(z) = / XB(N, a /L)(v~ 1 z)fi L (vN)di> 

Jueso(d+i) 

= —#(Z(L)nB(z,a/L)) 

and 

(o- * XB(N,a/L))(z) = a(B(z, a/ L)) 
and finally for any e > there is a big a such that 

(6) (1 - e)(or * XB(N,<x/L))(z) < (VL * XB(N,a/L))(z) < (1 + e)(flr * XB(JV,a/i))( 2; )' 

for any 2 and L > 1. 
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We take an arbitrary spherical cap B(z, r). We want to check that /j, L (B(z, r)) — > a(B(z, r)) 
as L — > oo. We fix an £ > and we take the convolution of (6) with the function ^vdt^tfv; 

v 7 cr(B(N,a/L)) 

with a very big L and this proves 

(1 - e)(flr * XB(J\T,r-a/X))(«) < (/iL * XB(N,r+a/L) ) (z) 

and 

Ol * XB(JV,r-a/i))(^) < (1 + e){a * XB(N,r+a/L))(z), 

for any z and L big. We take limits as L — > oo and since 2 is uniformly separated that means 
that [a l (B(z, r + a/L)\ B(z, r — ol/L)) — > uniformly in z G S d as L — > oo. Thus 

lim fi L (B(z,r)) = a(B(z,r)). 

L— »oo 

This is for an arbitrary spherical cap. This already implies the vague convergence of the mea- 
sures, see [Blu90], i.e. 



i y 



3=1 

for any / G C(§ d ). 

From the uniform densities condition on the Fekete points we may obtain other geometric 
consequences on the distribution of the Fekete points Z(L). We give one example. It is well 
known that using the bound given by by L. Fejes Toth [FT49] for the maximum of the minimal 
spherical distance between any set of tt l = (L + l) 2 points on § 2 , there exist at least two points 
in Z(L) at distance with 

j ^ cot 2 cj L -l (L + l) 2 vr 

d L < arccos , ujl ~ 



but as 



one has 



2 ' ^ (L + l) 2 -2 6' 
L arccos Ul f k = 3.80925 . . . , when L — > oo 



mmd(z LU z Lj ) < —. 
#3 L 



However the numerical results in [SW04] suggest that the right bound should be 

lim sup Lzmnd(zLi, zlj) < tt- 

To bound the maximal number, N, of disjoints spherical caps in S 2 of radius rj/L to be found 
in a larger spherical cap of radius a/Lwe use the following result due to J. Molnar 

N< * °(B(z,f)) 
~ VT2a(B(z,^)Y 
see [Mol52]. Then substituting in the density condition one gets 



77 < 4, /-^= = K. 

12 
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So, we get not only that there exists a couple of points in each generation Z(L) at distance 
smaller than k/L but a more uniform estimate: for any e > there is an a such that: 

z Li ,z Lj £B(z,a/L) Li 

for any spherical cap B(z, a/L). 

Up to now we have drawn information from the MZ- arrays and interpolating arrays to get new 
information on the Fekete points, but the reverse trend can also be useful. For instance, since 
any Fekete array has density one and a small perturbation makes it a MZ-array (or interpolation), 
then we obtain the following corollary 

Corollary 3.3. Given any e > 0, there are arrays Z £ and Z_ £ with densities D + (Z £ ) = 
D~(Z £ ) = 1 + e, D + (Z^ £ ) = D~(Z^ £ ) = 1 — e, such that Z £ is an L P -MZ array for any 
p £ [1, oo] and Z_ £ is an L p -interpolating array for any p £ [1, oo]. 

Thus the necessary density conditions that A/Z-arrays and interpolating arrays satisfy are 
sharp. 
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